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Effects  of  Numbers,  Sizes  and  Crystal 
Structure  on  Measurements  of  Ice  Hydrometeors 


1.  INTRODUCTION 


Measuring  numbers  and  sizes  of  Ice  crystals  has  always  been  a  difficult  and  frustrating 
task.  However,  the  accurate  description  of  Ice  hydrometeor  environments  continues  to  be 
Important  In  a  number  of  Investigations  such  as  those  concerned  with  electromagnetic 
attenuation  and  missile  erosion. 

The  early  days  of  Ice  crystal  measurements  Involved  long  tedious  hours  of  counting 
particle  Imprints  on  coated  glass  slides  and  film  or  on  aluminum  foil  that  had  been  exposed  to 
Ice  crystal  populations.  The  Introduction  of  the  PMS  Instruments*  In  the  early  1970's  greatly 
simplified  the  acquisition  of  hydrometeor  data  by  providing  the  means  to  electronically  count 
and  size  airborne  particles.  These  measurements  produce  number  versus  size  distributions 
that  describe  particular  hydrometeor  situations.  Thus,  the  problems  associated  with  Ice 
crystal  determinations  has  now  shifted  from  one  of  data  gathering  to  one  of  data 
»  Interpretation. 


(Received  for  publication  27  March  1990) 

*  Knollenberg,  R.G.  (1970)  The  optical  array:  an  alternative  to  scattering  or  extinction  for 
airborne  particle  size  determination,  J.  AppL  Meteorol  1:  86-103. 
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Typical  measurements  of  preclpltable  Ice  hydrometeors  by  a  PMS  one-dlmentlonal  probe 
(1-D)  show  the  largest  number  concentration  at  some  small  particle  size  with  the  numbers 
decreasing  as  the  sizes  of  the  hydrometeors  Increase.  When  the  numbers  are  plotted  versus  size 
In  a  seml-logarlthmlc  format,  that  decrease  Is  found  to  be  systematic  and  can  be  approximated 
by  a  straight  line.  Thus,  number  versus  size  distributions  may  be  described  by  exponential 
functions.  This  phenomenon  was  also  found  In  measurements  of  rain  and  was  first  reported 
by  Marshall  and  Palmer^  in  1948.  The  similarities  In  the  shapes  of  measured  number  density 
distributions,  both  liquid  and  Ice,  gives  rise  to  the  possibility  that  all  preclpltable 
hydrometeors  may  naturally  exist  In  exponential  form. 

There  have  been  however,  several  recurring  Inconsistencies  In  Ice  crystal  data  from  the  1-D 
precipitation  probe  that  are  difficult  to  explain,  especially  for  smaller-sized  hydrometeors. 

For  Instance,  although  the  major  portions  of  those  type  distributions  may  be  described  by 
exponential  functions,  the  first  or  first  few  class  sizes  frequently  show  deficient  number 
counts.^  One  possible  explanation  Is  that  evaporation  may  be  the  cause  of  the  lesser  number  of 
smaller  particles. 

Another  problem  Is  non-conforming  number  counts  In  the  larger-sized  classes  of 
precipitation  that  generally  Indicate  concentrations  larger  than  that  predicted  by  the 
exponential  trend  formed  by  the  majority  of  points. 

Recurring  disturbances  are  also  present  In  the  measurements  recorded  by  the  cloud  probe. 
Plots  of  those  number  densities  display  many  varied  shapes  and  often  have  values  In  their 
larger-sized  classes  that  do  not  conform  with  those  from  the  rest  of  the  distribution  and/or 
with  those  from  the  precipitation  measurements.  Also.  In  many  cases,  there  are  disparities  In 
the  normalized  numbers  (that  Is,  numbers  |>er  millimeter  bandwidth)  In  the  overlapping  size 
region  of  the  cloud  and  precipitation  Instruments. 

Figure  1  Is  a  number  density  plot  that  shows  examples  of  the  Inconsistencies  mentioned 
above.  The  data  were  recorded  on  15  February  1979  at  an  altitude  of  5562  meters  and 
temperature  of  -29.2'’C.  The  predominant  crystal  type  was  Identified  as  being  columns. 


^  Marshall,  J.S.  and  Palmer,  W.M.  (1948)  The  distribution  of  raindrops  with  size,  J.  Meteorol 
5:165-166. 

^  House,  R.A.,  Hobbs,  P.V.,  Herzegh,  P.H.,  and  Parsons,  D.B.  (1978)  Airborne  measurements  of 
the  size  distributions  of  precipitation  particles  In  frontal  clouds,  Conf.  on  Cloud  Physics  and 
Atmospheric  Electricity,  Jul  31-Aug  4.  Issaquah.  Wash. 
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Figure  1.  Processed  PMS  1-D  Data  Taken  on  15  Feb  1979  at  an  Altitude 
of  5562  Meters  and  Temperature  of  -29.2'’C. 


A  study,'*  conducted  In  1981,  ejqjlored  the  measuring  of  various  shapes  of  ice  ciystal  forms 
by  a  PMS  Instrument.  The  data  were  simulated  specifying  distributions  of  numbers  and 
sizes  and  allowing  each  individual  crystal  to  have  a  randomly- selected,  spatial  orientation.  It 
was  found  that  the  measurements  could  be  substantially  different  than  those  from  the  initial, 
actual  distributions  especially  in  cases  of  elongated  crystalline  structures. 

The  computer  program  used  to  fabricate  single  crystal  measurements  has  since  been 
expanded  to  allow  the  combining  of  individual  crystals.  The  simulated  data,  in  class  sizes 
common  to  those  of  the  PMS  Instruments,  resulted  in  distributions  similar  in  form  to  those 


^  Berthel,  R.O.  (1981)  The  Conversion  of  Aircraft  Ice  Crystal  Measurements  into  Terms  of  Liquid 
Water  Using  Simulated  Data.  AFGL-TR-8 1-0173,  ADA106419. 
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from  In-sltu  measurements.  Different  assumptions  of  crystal  numbers  and  sizes  were  found  to 
affect  those  distributions  In  varying  degrees  and  could  reproduce  the  Inconsistencies  cited 
above. 

This  report  demonstrates  the  effects  that  numbers  and  sizes  of  both  single  and 
combinations  of  randomly  oriented,  cylindrical-columnar  crystals  may  have  on  PMS  1-D 
measurements.  It  shows  that  measured  number  densities  can  vary  because  of  crystal  structure 
and  concentration  with  the  values.  In  some  cases,  having  little  resemblance  to  the  sizes  of  the 
particles  being  sampled.  This  study  Is  descriptive  In  nature  only  and  does  not  delve  Into  the 
conversion  of  Ice  hydrometeor  measurements  to  equivalents  of  liquid  water  that  are  necessary 
for  calculations  of  liquid-water  contents.  As  such.  It  Is  Just  the  first  step  In  understanding  the 
Inherent  complexities  Imbedded  In  measurements  of  Irregular-shaped  Ice  hydrometeors  and 
for  the  subsequent  Interpretation  of  those  measurements  In  the  realistic  definition  of  Ice 
crystal  environments. 


2.  SIMULATED  MEASUREMENTS 


The  data  used  In  this  Investigation  are  composed  of  simulated  measurements  of  columnar 
crystal  forms  of  known  numbers  and  sizes  with  PMS  1-D  probes. 

The  1  -D  instrument  (simply  stated)  consists  of  a  series  of  photodiodes  Illuminated  by  a 
laser.  A  hydrometeor  passing  through  the  sampling  volume  occludes  one  or  more  diodes 
thereby  producing  a  measurement  of  one  unspecified  dimension  that  Is  referred  to.  In  this 
rejx)rt,  as  a  shadow  length  or  L^.  The  output  from  the  aneilyses  of  the  1-D  data  consists  of  the 
numbers  of  hydrometeors  contained  In  15  size  classes  with  limits  predetermined  by 
Instrument  construction  and  calibration.  In  this  Investigation,  the  assumed  measurements  of 
the  cloud  probe  spans  the  range  from  0.01  to  0.31  mm  In  class  widths  of  0.02  mm.  The 
precipitation  probe  covers  0.2  to  4.7  mm  In  class  widths  of  0.3  mm. 


2.1  Single  Crystals 

The  assumed  ice  crystals  used  In  this  study  are  solid  columns  with  lengths  (L)  four  times 
their  diameters  (D).  If  one  were  to  look  directly  down  upon  such  a  falling  column  as 
diagrammed  In  Figure  2,  the  horizontal  plane  would  be  that  of  the  page  surface  delineated  by 
the  X  and  y  axes  where  the  vertical  plane  Is  perpendicular  to  that  surface,  defined  by  the  z  axis. 
The  diode  array  Is  assumed  to  be  aligned  with  and  below  the  paper's  surface  with  the  laser 
beam  emanating  from  directly  above.  A  column,  with  L=1  mm,  passing  through  a  1-D 
Instrument  with  Its  longitudinal  axis  aligned  with  the  diode  array  In  both  the  horizontal  and 
vertical  planes,  as  In  Figure  2a,  would  result  In  a  L,  measurement  of  the  true  length  or  1  mm. 

If  however,  the  colurrm  were  rotated  90°  in  the  horizontal  plane  as  In  Figure  2b,  the  would 
be  equal  to  D  or  0.25  mm. 
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a.  Column  aligned  with  b.  Column  of  la  c.  Column  of  la 

diode  array  In  both  rotated  90°  In  rotated  45°  In 

horizontal  and  vertical.  horizontal  plane.  vertical  plane. 


Figure  2.  Diagrammatic  Representation  of  the  Measuring  of  a  Column 
with  a  Diameter  of  0.25  mm  and  Length  of  1  mm  In  Three  Particular 
Spatial  Orientations. 


The  Lg  Is  further  affected  by  the  crystal's  orientation  In  the  vertical  plane  as  any  degree  of 
tilt  of  the  column's  longitudinal  axis  from  that  of  the  diode  array  In  the  z  direction  will  result 
in  a  smaller  shadow  length  as  demonstrated  In  Figure  2c.  The  1  mm  column  of  Figure  2a, 
aligned  with  the  diode  array  In  the  horizontal  plane  but  now  at  an  angle  of  45°  In  the  vertical, 
will  give  a  Lg  of  0.707  mm.  Thus,  a  crystal's  spatial  orientation  In  both  the  horizontal  and 
vertical  planes  has  a  direct  effect  on  the  measurement  obtained  with  a  1-D  type  Instrument. 

This  Investigation  assumes  a  hypothetical  scenario  where  PMS  probes  are  being  flown 
through  a  known  Ice  hydrometeor  environment  and  are  measuring  a  specific  number  of 
columns  of  defined  sizes.  Thus,  the  simulation  of  data  begins  with  the  sijeclficatlon  of 
numbers  and  sizes  of  columnar  crystals.  Each  column  Is  then  assigned  a  randomly  selected 
angle  to  define  the  position  of  Its  longitudinal  axis  with  respect  to  a  probe's  diode  array. 

Ninety  degrees  of  freedom  are  allowed  In  the  horizontal  plane  and  45°  In  the  vertical.  Forty- 
five  degrees  of  freedom  In  that  plane  are  considered  reasonable  for  environments  of  average 
wind-shear  turbulence  whereas  90°  would  represent  severe  conditions  and  tumbling  crystals. 
The  resulting  Lg  Is  then  placed  In  the  appropriate  class-size  category  thereby  mimicking  the 
measuring  of  a  crystal  by  a  1-D  Instrument. 

Again,  It  must  be  emphasized  that  this  Investigation  Is  Intended  to  be  descriptive  only.  No 
consideration  has  been  given  to  Instrument  calibration,  sampling  times  or  aircraft  speed. 
Therefore,  the  numbers  cited  in  the  simulated  measurements  will  be  In  general  terms  (numbers 
per  unit  volume  per  probe  bandwidth,  N  V“* 
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Figure  3  shows  the  results  of  the  fabricated  measurements  of  1000  single  columns  with 
L=1  mm  In  the  class  widths  previously  specified  for  the  assumed  1-D  probes  and  In  a  seml- 
logarlthmlc  format  common  for  number-density  distribution  plots.  If  each  of  these  columns 
were  precisely  aligned  with  the  1-D's  diode  array,  all  of  the  1000  measurements  would  have 
been  contained  In  the  third  precipitation  class  (0.8  to  1.1  mm).  Because  of  the  crystal  structu^"'? 
and  random  spatlail  position,  Just  ~45  percent  were  recorded  as  being  In  that  size  range  with 
the  remainder  being  measured  as  smaller  particles.  Notice  that  some  of  the  crystals  were  i 

oriented  In  such  a  manner  that  the  values  could  also  be  measured  by  a  cloud  probe.  The 
numbers  and  sizes  shown  In  Figure  3.  and  In  all  the  simulated  measurements  In  this  report, 
represent  typical  examples  only,  since  random  positioning  will  produce  variations  In  separate 
simulations  of  Identical  crystals. 
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Figure  3.  Simulated  Measurements  of  1000  Columns  with  Diameters 
of  0.25  mm  and  Lengths  of  1  mm. 


6 


2.2  Combinations  of  Crystals 


When  the  single  columns  described  above  are  allowed  to  attach  to  one  another  In  a  random 
fashion,  the  resulting  combinations  display  varying  shapes  from  elongated  chains  to  rosettes 
or  tightly  packed  clusters.  Figure  4  shows  computer  constructed  combinations  of  2,  3,  4,  and  5 
columns  In  several  different  configurations.  Viewing  these  diagrams  makes  It  easy  to  envision 
the  ntyiiad  of  measured  sizes  that  may  result  from  differences  In  crystal  spatial  orientation 
and  position  of  attachment. 


Figure  4.  Computer  Generated  Examples  of  Combinations  of  2,  3,  4. 
and  5  Columns. 
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Figure  5  shows  the  simulated  measurements  of  1000  of  each  of  these  combinations  In  the 
same  seml-logaiithmlc  plotting  format  as  the  single  columns  In  Figure  3.  The  random 
selection  of  column  positions  and  attachments  tends  to  produce  forms  that  approximate 
normal  curves  with  both  the  peak  number  values  and  range  of  sizes  Increasing  as  the 
combinations  become  larger.  Also,  the  numbers  In  the  larger  sized  classes  form  general 
exponential  shapes.  There  are  no  counts  attributed  to  the  cloud  probe  tn  these  four  cases 
although  the  possibility  does  exist  for  any  one  of  these  combinations  to  attain  the  necessary 
configuration  and  spatial  orientation  to  fit  into  one  of  the  last  three  classes  of  that  assumed 
instrument. 
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Figure  5.  Simulated  Measurements  of  1000  of  Each  Combination  of 
2,  3,  4,  and  5  Columns  with  1  mm  Lengths. 
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Five  columns  per  combination  are  the  maximum  considered  In  this  study  as  the 
mechanism  of  random  adhesion  used  In  this  model  tends  to  favor  the  stacking  of  crystals  auid 
produces  "snowball"  type  structures  from  combinations  consisting  of  larger  numbers. 


2.3  Mixtures  of  Single  Columns  and  Combinations  of  Columns 

j 


Although  measurements  of  small  sized  hydrometeors  are  sometimes  conducted  In 
environments  consisting  of  all  single  colutims  or  nearly  Identical  combinations.  It  Is  more 
common  to  find  singles  coexisting  with  various  combinations. 

Specific  populations  of  columns  containing  both  single  crystals  and  combinations  may  be 
defined  by  assigning  numerical  values  to  each  hydrometeor  category.  This  creates  an  assumed 
distribution  that  may  be  described  mathematically  where  the  concentration  Is  a  function  of 
the  number  of  columns  making  up  a  hydrometeor.  However,  It  Is  desirable  to  describe  the 
distribution  In  conventional  terms  of  numbers  versus  sizes,  except  that  the  sizes  of 
combinations  cannot  be  specified  because  of  uncertainties  associated  with  the  random 
attachment  of  crystals  In  the  building  of  the  hydrometeors.  To  overcome  this  ambiguity,  the 
numbers  of  assumed  hydrometeors  are  defined  as  exponential  functions  of  their  equivalent- 
melted  diameters  (Dg)  or.  In  other  words,  the  diameters  of  the  spherical  water  drops  resulting 
from  the  melting  of  the  column  or  columns  making  up  the  hydrometeor.  (The  numbers  and 
sizes  of  the  assumed  particles  subjected  to  simulated  measurement  will  be  designated  by  circles 
In  all  the  following  plots  except  In  those  cases  where  they  are  specified  In  the  figure  captions.) 

When  the  numbers  and  sizes  used  In  developing  Figures  3  and  5  (1000  single  columns  and 
1000  of  each  combination  of  2,  3,  4,  emd  5  columns  with  L=1  mm)  are  assumed  to  exist  as  a 
mixture,  the  assumed  distribution  of  hydrometeors  forms  an  exponential  slope  (A)  of  zero.  The 
simulated  measurements  of  that  mixture  are  shown  plotted  In  Figure  6.  The  values  in  the 
precipitation  classes  display  a  general  exponential  shape  for  the  latter  portion  of  the 
distribution  (as  per  the  solid  line  In  the  figure)  with  the  first  four  showing  a  deficiency  In 
numbers  with  respect  to  the  extrapolated  exponential  (dashed)  line.  The  deviation  from 
exponential  shape  Is  similar  to  but  more  pronounced  than  that  normally  observed  In  actual 
PMS  measurements  of  Ice  crystals.  The  plot  also  shows  a  small  number  of  counts  In  the  latter 
three  cloud  probe  classes  similar  to  that  of  Figure  3. 
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Figure  6.  Simulated  Measurements  of  an  Assumed  Distribution  of 
1000  Single  and  1000  of  Each  Combination  of  2,  3,  4,  and  5  Columns 
with  1  mm  Lengths. 


3.  EFFECTS  OF  CRYSTAL  SIZE 


When  the  L  of  the  Individual  columns  used  In  Figure  6  Is  changed  to  0.5  mm  and  the  same 
number  of  single  and  combinations  of  crystals  are  assumed  to  be  measured,  the  plot  of 
Figure  7  results.  In  this  case,  numbers  are  found  In  only  five  classes  of  precipitation  with  the 
last  four  being  In  exponential  form.  Although  a  difference  still  exists  between  the  number  In 
the  first  class  and  the  extrapolated  exponential  line,  the  deficit  Is  considerably  less.  The  slope 
of  the  exponential  line  over  the  larger  class  sizes  Is  also  much  steep>er. 
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Particles  able  to  be  measured  by 


Figure  7.  Simulated  Measurements  of  an  Assumed  Distribution  of  1000 
Single  and  1000  of  Each  Combination  of  2.  3,  4,  and  5  Columns  with 
0.5  mm  Lengths. 


Another  outstanding  feature  Is  that  the  values  associated  with  the  cloud  probe  are 
considerably  greater  with  counts  In  the  latter  10  classes.  The  numbers  c  ntalned  in  the  first 
four  of  these  classes  (108)  represent  particles  not  Included  In  the  precipitation  counts. 

Simulated  measurements  from  hydrometeors  of  columns  with  L's  between  0.5  to  1  mm 
show  that  the  numbers  In  the  smaller-sized  classes  of  the  precipitation  probe  deviate  further 
from  the  exponential  lines  as  the  columns  become  larger.  Also,  the  slopes  of  the  resulting 
distributions  become  more  shallow  with  larger  particle  size. 
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4.  EFFECTS  OF  RELATIVE-SIZE  NUMBER  CONCENTRATIONS 


Variations  In  measurements  caused  by  the  relative  numbers  of  different  sized  particles 
were  explored  using  two  scenarios.  First,  the  number  of  assumed  single  columns  was  large 
with  decendlng  amounts  as  the  crystals  combined  (negative  exponential  slope).  Second,  the 
number  of  singles  was  small  with  ascending  amounts  as  the  combinations  grew  larger 
(positive  exponential  slope). 

Figure  8  shows  the  simulated  1-D  measurements  from  an  assumed  hydrometeor 
distribution  of  single  and  combinations  of  1  mm  columns  where  numbers  vs  equivalent- 
melted  diameters  form  a  negative  exponential  slope  of  -10  mm“*.  The  last  five  classes  follow  a 
general  exponential  shape  with  the  first  two  showing  an  abrupt,  negative  departure  from  the 
stralght-Ilne  projection.  The  small  number  of  counts  In  the  last  three  classes  of  the  cloud 
probe  are  again  present  as  In  Figure  6. 
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Figure  8.  Simulated  Measurements  of  an  Assumed  Distribution  of  Single 
and  Combinations  of  1  mm  Columns  Whose  Numbers  Versus  Ekjulvalent- 
Melted  Diameters  Form  a  -10  mm“*  Exponential  Slope. 
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The  plot  of  Figure  9  also  exhibits  an  exponential  trend  In  the  latter  four  precipitation 
classes  but  over  larger  sizes.  The  Initial  assumed  distribution.  In  this  case,  was  defined  by 
a  positive  slope  of  10  The  first  few  classes  however,  display  an  exaggerated  departure 

from  the  projected  exponential  line  when  compared  to  real  data.  Also,  there  are  fewer  counts 
in  only  two  classes  of  the  cloud  Instrument. 
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Figure  9.  Simulated  Measurements  of  an  Assumed  Distribution  of  Single 
and  Combinations  of  1  mm  Columns  Whose  Numbers  Versus  Equivalent- 
Melted  Diameters  Form  a  10  mm  Exponential  Slope. 


It  Is  Interesting  to  note  the  similarities  In  the  range  of  class  sizes  that  form  the 
exponential  lines  In  the  plot  of  Figure  6.  which  represents  an  Initial  distribution  with  zero 
slope,  and  Figure  9  with  a  slope  of  10  mm”>.  The  plot  of  Figure  8,  with  an  Initial  slope  of 
-10  mm~',  forms  a  similar  exponential  but  over  a  range  of  smaller  sizes.  It  also  shows  more 
realistic  numbers  In  the  small  classes. 
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Simulated  measurement  were  also  made  on  a  relative- size,  hydrometeor  distribution  of 
A=-10  mm“*  composed  of  columns  with  L=0.5  mm.  The  results  are  shown  In  Figure  10.  When 
compared  with  Figure  7  where  A=0  mm“*,  this  plot  shows  the  counts  In  the  first  precipitation 
class  to  be  closer  to  the  exponential  projection.  In  fact,  the  counts  In  the  first  class  of  Figure  7 
are  smaller  than  that  of  the  second,  whereas  the  opposite  Is  shown  In  Figure  10.  The  numbers 
In  the  latter  three  classes  form  an  ejqxtnentlal  nearly  Identical  to  the  one  In  Figure  7  but  over 
a  smaller  range  of  sizes.  The  cloud-sized  counts  are  similar  to  Figure  7. 
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Figure  10.  Simulated  Measurements  of  an  Assumed  Distribution  of  Single 
and  Combinations  of  0.5  mm  Columns  Whose  Number  Versus  Ekjulvalent- 
Melted  Diameters  Form  a  -10  mm"^  Exponential  Slop)e. 


The  distributions  In  Figures  6  through  10  indicate  that  the  number  deviations  In  the  first 
and  first  few  classes  of  the  simulated  precipitation  data  are  def)endent  upon  both  the  sizes  and 
relative  numbers  of  the  small  hydrometeors.  The  resulting  slopes  of  the  latter  precipitation 
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classes,  on  the  other  hand,  seem  only  to  be  governed  by  larger  hydrometeor  sizes.  Thus,  any 
condition  that  affects  the  sizes  will  cause  a  variance  In  slope.  For  example,  adding 
combinations  consisting  of  larger  numbers  of  columns,  hence  larger  hydrometeor  sizes,  will 
give  lesser  slop)es  while  subtracting  causes  the  opposite  effect. 


5.  EFFECTS  FROM  A  MIXTURE  OF  CRYSTAL  SIZES 


If  a  hydrometeor  environment  were  composed  of  building  blocks  of  different  sized  Ice 
crystals,  one  would  expect  the  measurements  to  display  a  mixture  of  the  predominant 
characteristics  of  both  the  largest  and  smallest  particles.  Figure  11  shows  the  simulated 
measurements  from  a  rough  approximation  of  such  a  situation. 

The  initial  distribution  In  this  example  was  formulated  from  0.5,  0.75,  and  1  mm  columns 
In  singles  and  combinations  as  In  the  previous  plots.  The  number  versus  slope  was 
assumed  to  be  -10  mm'*  thus,  the  results  of  the  simulated  measurements  may  be  compared  to 
those  plots  that  correspond  to  the  largest  (Figure  8)  and  smallest  (Figure  10)  column  sizes. 

In  this  case,  the  first  precipitation  class  has  slightly  more  counts  than  the  second,  very 
similar  to  Figure  10  but  has  a  much  smaller  number  discrepancy  from  the  projected 
exponential  line.  This  Is  because  of  the  more  shallow  slope,  approximately  that  of  Figure  8. 
As  expected,  the  number  of  particles  with  Lg  sizes  that  could  be  measured  with  a  cloud  probe 
are  similar  to  those  of  Figure  10. 
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O  Numbers  and  D,  sizes  of 
assumed  hydromefeors 


Figure  1 1 .  Simulated  Measurements  of  an  Assumed  Distribution  of  Single 
and  Combinations  of  Columns  with  0.5,  0.75  and  1  mm  Lengths  Whose 
Number  Versus  Equivalent-Melted  Diameters  Form  a  -10  mm"*  Exponential 
Slope. 


6.  EFFECTS  FROM  SIMULTANEOUS  HYDROMETEOR  MEASUREMENTS 


The  inconsistencies  normally  present  in  the  largest-sized  measurements  of  the  1-D 
precipitation  probe  are  apparently  caused  by  a  relatively  small  number  of  hydrometeors. 
Because  the  seml-logarithmlc  number  versus  size  distribution  plot  has  small  counts  tn  the 
larger-sized  classes,  the  addition  of  one  or  two  particles  may  cause  a  wide  disparity  from 
exponential  shape.  The  random  orientation  of  elongated  ice  hydrometeors  may  be  an 
explanation,  although  this  study  indicates  that  a  very  low  probability  exists  for  the  formation 
of  such  structures  with  the  necessary  spatial  positioning.  However,  if  it  is  possible  for  two 
h3rdrometeors  to  form  a  loose  attachment  and  pass  through  the  1-D  instrument 
simultaneously,  then  this  phenomenon  would  be  the  anticipated  result. 
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A  simulation  of  this  effect  Is  shown  In  Figure  12.  The  X's  In  this  plot  are  the  simulated 
precipitation  probe  measurements  from  columns  of  L=0.7  mm  with  single  crystals  and 
combinations  consistent  with  the  previous  plots.  The  Initial  distribution  has  a  slope  of 
-10  mm”'.  The  line  shows  the  exponential  trend  formed  by  the  last  four  classes. 


9 

1 

z 

05 

z 

a 

oc 

Ui 

a 

lit 

2 
D 


Z 

D 

cc 

lit 

a. 

oc 

lit 

CQ 

2 

O 


Particles  able  to  be  measured  by 


MILLIMETERS 


Figure  12.  Simulations  Showing  the  Possible  Effects  from 
Simultaneous  Measurements  of  Hydrometeors. 


Simultaneous  measurements  were  simulated  by  randomly  selecting  20  values  from  this 
distribution,  subtracting  them  from  their  classes  and  forming  them  Into  10  pairs.  The 
constituents  of  each  pair  were  then  assumed  to  be  attached  with  new  sizes  equal  to  the  sums  of 
the  Individual  L,  values.  Classifying  these  combined  particles  thus  produced  an  adjusted 
distribution  with  a  net  change  in  total  number  of  just  10  hydrometeors  or  -0.5  percent. 
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The  results  of  three  such  adjustments  displaying  noticeable  changes  In  the  larger 
precipitation  class  sizes  are  plotted  In  Figure  12  and  are  designated  by  circles,  squares  amd 
crosses.  Distribution  shapes  similar  to  these,  where  the  values  of  the  larger  precipitation 
classes  deviate  from  the  exponential,  are  not  uncommon  In  PMS  1-D  measurements  (see 
Figure  1). 


7.  EFFECTS  ON  THE  OVERLAPPING  REGION  OF  THE  CLOUD  AND 
PRECIPITATION  PROBES 


The  physical  nature  and  random  positioning  of  Ice  crystals  may,  as  demonstrated  for 
preclpltable  particulates  In  the  preceding  sections,  affect  the  measurements  from  both  the  PMS 
cloud  and  precipitation  probes.  The  severity  of  the  effect  Is  associated  with  hydrometeor 
shap)es,  and  sizes  of  the  particles  and  class  widths  of  the  measuring  Instrument.  These 
phenomena  may  help  to  explain  some  of  the  Inconsistencies  observed  In  the  overlapping  size 
region  of  the  cloud  and  precipitation  Instruments. 

Figures  6-11  show  that  elongated  particles,  because  of  spatial  positioning,  may  be  placed  In 
smaller  class  sizes  than  their  lengths  would  Indicate.  Thus,  a  particular  particle  with  a  length 
normally  associated  with  the  precipitation  probe  may  actually  be  recorded  by  the  cloud 
Instrument.  Measurements  from  the  cloud  probe  would  also  be  affected  in  the  same  manner. 

In  addition,  although  both  probes  would  be  exposed  to  the  same  hydrometeor  environment 
during  actual  measurements,  they  would  be  sampling  different  particles  of  various  shapes  In 
different  spatial  positions.  This  adds  another  degree  of  uncertainty. 

If  measurements  were  made  In  an  environment  that  contained  Ice  spheres  with  sizes 
spanning  both  those  of  the  cloud  and  precipitation  probes,  the  measurements  from  the  two 
probes,  when  "normalized"  (divided  by  Instrument  bandwidth  to  give  numbers  per  millimeter 
bandwidth),  should  form  a  smooth  transition  with  number  densities  that  accurately  describe 
the  hydrometeor  population.  If  the  Ice  spheres  were  changed  to  columns,  the  normalized 
values  may  be  compatible  but  will  not  be  true  measurements  of  the  hydrometeor  numbers  and 
sizes.  Because  of  the  uncertainties  associated  with  measurements  of  different  parcels  of  air 
with  separate  Instruments  and  from  random  crystal  attachment  and  spatial  positioning, 
perfect  compatibility  should  not  be  ejqiected.  This  type  situation  would  be  found  m  generating 
cloud  systems  where  Ice  crystals  are  being  formed. 

In  non-generating  cloud  systems,  smooth  transitions  In  the  measurements  from  the  cloud 
and  precipitation  probes  would  be  rare.  The  number  density  distributions  from  the  cloud 
probe  measurements  may  exhibit  many  varied  shapes  depending  upon  the  Individual 
populations  of  the  cloud  and  unrelated  preclpltable  hydrometeors. 

An  example  of  one  possible  effect  Is  shown  In  Figure  13.  In  this  scenario,  the  hydrometers 
producing  the  simulated  measurements  shown  In  Figure  10  are  assumed  to  be  falling  through 
non-generating  clouds.  The  Xs  are  the  normalized  values  from  the  precipitation  probe.  The 
normalized  cloud  probe  values  attributed  to  the  assumed  cloud  hydrometeors  are  designated  by 
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the  triangles.  If  the  numbers  of  preclpltable  particles  that  were  oriented  such  that  they  could 
be  measured  by  a  cloud  probe  were  added  to  the  assumed  cloud  hydrometeor  values,  the 
resulting  normalized  distribution  would  be  that  shown  by  the  crosses.  Distributions  from  In- 
sltu  measurements  frequently  show  similar  shapes  (see  Figure  1). 
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Figure  13.  Simulated  PMS  1-D  Measurements  of  Assumed  Cloud 
Hydrometeors  In  Concert  with  the  Preclpltable  Particles  Assumed 
In  Figure  10.  This  hypothetical  situation  shows  an  example  of 
Incompatibility  between  cloud  and  precipitation  probe  data. 
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8.  CONCLUSIONS 


The  Investigation  described  In  this  report  used  simulated  measurements  of  assumed  Ice 
crystal  environments  to  explore  several  Inconsistencies  that  are  often  present  In  PMS  1-D 
recordings.  The  results  strongly  Indicate  that  the  spatial  orientation  of  Irregular  or  elongated 
particles  may  be  the  cause  as  they  can  affect  both  number  and  size  measurements. 

Several  Interesting  observations  surfaced  during  the  course  of  this  study.  For  Instance,  the 
simulated  data  show  that  measuring  and  classifying  randomly-oriented.  Irregular  particulates 
can  produce  an  apparent  deficiency  of  numbers  In  the  first  or  first  few  classes  of  the 
precipitation  measurements.  The  magnitudes  of  those  deficiencies  are  related  to  ciystalllne 
structure,  particle  sizes  and  numbers  of  those  hydrometeors. 

The  major  portions  of  the  simulated  precipitation  measurements  also  form  exponential 
shapes,  as  do  actual  data,  where  the  slopes  seem  to  depend  on  the  numbers  and  sizes  of  the 
larger  particles.  The  simulations  of  assumed  Ice  distributions  of  Identical  sizes  but  Avlth 
distinctly  different  number  concentrations  resulted  In  plots  showing  similar  exponential 
slopes  but  over  different  size  ranges.  However,  the  assumption  of  more  large  hydrometeors 
than  small  produced  a  distribution  shape  that  showed  an  exaggerated  departure  from  the 
exponential  line  In  the  small-sized  classes  when  compared  to  actual  measurements.  This 
would  seem  to  Indicate  that  natural  hydrometeor  populations  contain  larger  amounts  of 
smaller  particles  and  may,  in  fact,  form  exponential  distributions. 

This  study  also  showed  that  the  number  counts  in  the  cloud  probe  can  be  adversely  affected 
by  the  spatial  orientation  of  Ice  crystals.  Thus,  numerous  variations  are  p)ossible  In  both  the 
cloud  and  precipitation  measurements  depending  upon  the  particular  environment  at  the  time 
of  sampling.  On  the  other  hand,  the  Inconsistent  scattered  number  counts  In  the  larger  class 
sizes  of  the  precipitation  probe  may  be  caused  by  the  combining  of  hydrometeors  or  by  the 
simultaneous  measuring  of  more  than  one  particle. 

The  columnar  structure  assumed  In  this  investigation  was  chosen  because  It  Is  often 
Identified  as  the  prevalent  form  In  small  ice  crystal  environments.  The  random  combining  of 
these  typje  crystals  also  tend  to  form  rosettes  or  combinations  similar  to  those  found  In 
nature.® 
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